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Abstract 


Xn this paper we Introduce the notion oi radio Jets turbulent mixing 
regions. We further propose that the essential small scale viscous dissipa- 
tion In these Jets Is by emission of MHD waves and by their subsequent strong 
damping due, at least partly, to gyro-resonant acceleration of supra-thermal 
particles. The equilibrium eddy, wave and particle spectra are not found 
exactly In this paper, but the problem Is defined and rough estimates of the 
spectra are given to aid In the observational Interpretation. 

A formula relating the synchrotron surface brightness of a radio Jet to 
the turbulent power Input Is deduced on very general and simple grounds from 
our physical postulates, and it Is tested against the data for NGC315 and 
3G31 (NGC383). The remarkable predictive power of this formula (Section III 
of this paper) Is the major Justification for presenting our model at this 
time. The predicted brightness depends essentially on the collimatlon 
behavior of the Jet, and, to a lesser extent, on the CH picture of a "high" 
nozzle with accelerating flow. The conditions for forming a large scale Jet 
at a high nozzle from a much smaller scale Jet are discussed. In particular 
the optimum condition for retaining the memory of the initial Jet direction 
Is given as that of a turbulent Jet. 

The effect of entrainment on our prediction Is discussed with the use of 
similarity solutions. Although entrainment Is Inevltrbly associated with the 
turbulent jet. It may or may not be a dominant factor depending on the ambient 
density profile. 
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I. INTRODUCTION 


It is becoming dear that relativistic particlu reacceleration occurs 
at distances £rom 1 to 40 kpc from the active nuclei of some radio galaxies, 
in the rapidly-expanding bases of radio jets linking these nuclei to their 
distant radio lobes. The evidence for this particle reacceleration is of 
three kliids. First, the observed surface brightnesses of radio Jets 
decrease with Jet width much less rapidly than expected if they were adlaba- 
tically-expanding fluxes of particles and fields without particle replenish- 
ment or field amplification. This effect has been found in most of the well- 
studied radio jets (3C31: Burch 1979, Fomalont et al. 1980; NGC315: 

Willis ^ al^ 1981a; HB13; Masson 1979; 3C296 and 3G449: Birkinshaw ®t al. 
(1981). Second, a polarized optical continuum coincides with the bright radio 
knots in the Jets of M87 (Owen et al. 1980) and 3C277.3 (Bridle et al. 1981). 
If this optical continuum is Indeed synchrotron radiation, the lifetimes of 
the radiating electrons in the equlpartltlon fields in these knots are much 
less than the light travel times to the knots from the galactic nuclei. This 
is clear evidence for relativistic particle replenishment within the knots. 
Third, the longer radio Jets show little or no gradient in radio spectral 
index along their lengths (e.g. Burch 1979; Willis et al. 1981 a,b). The 
evidence for large-scale particle reacceleration from such radio spectral 
data is more dependent on assumptions about particle transport velocities in 
the Jets than the first two lines of evidence, but is nevertheless suggestive. 

In this paper we report excellent fits to the observed synchrotron 
brightness distributions of the well-observed radio Jets in 3C31 and NGC315, 
based on a simple theory derived from the scaling laws applicable to "fully 
developed" volume turbulence. This turbulence is anisotropic because the 
mean velocity of a Jet imposes a preferred direction on the flow, but we 
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suppose that it may be treated as homogeneous and isotropic in the co-movlng 
frame of the Jet. 

The volume turbulence invoked in our model is generated in the free mixing 
layer that develops following the emergence of a jet with large Reynolds 
number from a high "nozzle" (e.g. chan and Henriksen 1980 (CH) ; Bridle# Chan 
and Henriksen 1981 (BCH)) or when the jet encounters a density "plateau" in 
its path (see Section V below) . The essential features of this turbulence for 
our purposes are that it exhibits large scale structure (e.g. Brovm and Roshko 
1974? Brown and Thomas 1979; Cantwell 1981, for a review) and soon fills the 
entire voluine of the jet, eliminating the laminar (potential) core altogethvsr. 
The growth of this large scale structure may be understood at each stage 
in terms of ring vortex pairing, mutual orbiting and merging followed by 
these processes repeating with the Just-merged eddies on a larger scale. 

(This is really a description of the non-linear growth of a Kelvin-Helmholtz 
instability). During this development, ambient material is entrained and 
Intense smaller scale turbulence is generated in the regions between the 
vortices (Cantwell 1981) , presumably establishing the turbulent cascade to 
higher wavenumbers which is eventually dominated by viscosity on the Kolmo- 
gorov microscale, . The small scale of the stretched turbulence generat- 
ing regions between the vortices can be associated with the Taylor microscale, 
(Cantwell 1981). 

We envisage that a primarily hydrodynamic turbulent cascade exists in a 
radio jet (a) because the Debye length is small (?: 2xl0®cm) and yet there are 
many particles in a Debye sphere (i 4x10^®), (b) because we expect the magnetic 
field to be dynamically weak (CH; BCH) except in the smaller eddies (CH? De 

Young 1980) and (c) because the small gyroradii of the particles ('v-O.l AU for 

-5 

fields 'v 10 gauss) should impart a fluid character to otherwise collisionless 
systems. The central idea which we use throughout for coupling the eddy 
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cascade to the MHD wave spectrum is the Llghthill MtID radiation from the 
turbulent eddies (e.g. Kulsrud 1955; Kato 1966; Stein 1981 for a summary). 

We further suppose that the viscous sub-range of the eddy cascade is estab- 
lished by this wave omission and subsequent damping, principally by particle 
acceleration. We defer to a later paper the detailed calculation of self- 
consistent eddy, wave and particle spectra that is required to Justify our 
vieW), proceeding here instead by plausibility arguments and by direct appeal 
to the observations of radio Jots (Section III of this paper). 

The radiated mhd waves may accelerate electrons directly by a resonant 
interaction (e.g. Laoomba 1977) or by developing into collisionless shock 
waves with scattering wakes which can mediate the first order Fermi acceleration 
process (e.g. Bell 1978a, b; Blandford and Ostriker 1978; Bicknell and Melrose 
1981) . Our estimates suggest that non-linear waves and shock waves can occur 
at I while true dissipation of the turbulent energy occurs in a range of 

smaller scales which are in resonance with typical high energy electrons 
present in Che Jets. 

This non-linear scheme for particle acceleration avoids the difficulty of 
propagating surface-generated resonant waves away from the limb of a jet 
encountered (Eilek 1981) by the pioneering idea of Pacholczyk and Scott (1976) . 
In our view, the strong damping of the resonant MHD waves is not a problem, but 
rather a statement of the strong coupling between particle energy and turbulent 
energy spectra to be expected locally . Section III below will give 
observational evidence for this view. Our model should also be distinguished 
from that of a turbulent wave cascade driven at the top end by Kelvin-Helm- 
holtz instabilities (see Benford, Ferrari and Trussoni, .1980; foe a recent 
view) . 

Section II of this paper derives (from simple scaling arguments) our basic 
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expression for the variation of the brightness B of a jet with jet radius R. 
Section III tests this relationship against observations of NGC 315 and 3C 31. 
Sections IV and V discuss entrainment and the relation between pc-scale and 
kpc-scale radio jets in the light of our general principles. Section VZ 
contains our conclusions and inoicates what remains to be done. 
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Hi RADIO BRIGHTNESS PROFILE OF TURBULENT JETS 
(a) Principles 

We aBsutoe lor dellnlteness that the apectrum of the turbulent inertial 
range is Kolmogorov (1968) In the frame co-movlng with the jet. It should be 
cleui: how to proce^id with a slight variation on the theme ^ however, such as a 
Kraichnan (MHD) spectrum (e.g. DeYoung 1980 and references therein). 

There is thus local homogeneity in the co-moving frame as well as a kind of 
transverse incompressibility (CH)i The range is taken to extend from 
I 'V R(Z) (Brown end Roshko 1974), the "radius” of the jet at a distance 
Z from a "nozzle" (or from a point where strong turbulent interaction with 
the surroundings begins), down to the effective damping or Kolmogorov scale 
. We suppose further that the "internal" or true dissipative viscosity 
in the jet is provided largely by the resonant interaction between IJHD waves 
ana the radiating electrons. This should be true even for the shock wave 
Fermi mechanism because the energy gained by the particles must derive 
ultimately from the turbulence through the mechanism of the pre-shock and 
post-shock tcatterers, of which the resonant waves are the most efficient. 

We note that we expect the turbulent "viscosity" (i.e. the Reynolds 
stress) to provide the effective "external" or entrainment viscosity, but 
this is not necessarily dissipative and is t, of course, to be used to 
calculate the Reynolds number. 

The Alfven wavelength in resonance with electrons may be written simply 


as (2 itw) r 


eg 


when the Alfven speed is much less than that of light and 

'mi 


the wave direction is close (i.e. « -~) to that of the mean magnetic 

r" 

field (e.g. Zilek, 1979). Here r is the mean electron gyro-radius and 

eg 

p is the cosine of the mean electron pitch angle. Hence the eddy scale 
which will couple to the electrons by the emission of these waves is 
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^ - f aiiM) "• 2TtY,m^cV(®i‘) 


( 1 ) 


where the factor f (<l) aUowa for the difference between the eeaXe of an eddy 
and ita characteristic Llghthill wavelength («ee below) » 5 is the spatially 
averaged niagnetlc field strength, and the other symbols have their uaual mean- 
ings. (Note that the effective damping scale is not uniquely defined in the 
absence of a self-consiatent description of eddy, wave and particle spectra. 

See also below.) Numerically, Eq. (1) gives 'V 10 ^‘*(y^( 5)/B(-5)) fp cm 
where the numbers in parentheses give the logarithmic unit of the stJmdard 

measure in e.g.s. units. The necessity for this latter scale to be considerably 

1 !%. 

greater than the ion gyroradlus ('vio^ Tg cm) but substantially less than the 
Taylor microscale suggests that the mechanism discussed here can accelerate 
electrons in the range 10^ < • We therefore use y ■ 10^ fo>^ our 

estimates. Eilek (1981) has shown that Alfven waves are absorbed on scales 
(10^ - lO**)! (the coiiffioient increases with X) for conditions appropriate 

to the radio jets we seek to model, so that S 10^ is also set by requiring 
that the waves are absorbed in less than a jet width ('vlO^^ cm) . 

The Reynolds number Re of the jet flow is fixed by assumption (1) to 
be (e.g. Landau and Lifshitz 1959) : 

Re » 2 x1q 9 (B(-5)R(21)/y (5))'*^^ (fp)"'*'^^ (2) 

Jv 6 

such a Reynolds number should be compared with those of 'vlO® obtaining in 
various shipwake problems (Cantwell 1981) ; our postulated range in scales 
(the jet radius R to £ ) is evidently not far removed from experience. 

One of the "a priori" arguments in favor of the observed radio jets 
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b«lng fcuiebulant is th« similarity of the obsorvod "con« angles" near the bases 
of various jets. For example, in each of NC5C 315 (WIHIb et al . 19f)la) , 3C 
31 (Bridle et al. 1980) and NGC 62S1 (Willis et al . 1981b) , dMIZ is initially 
'vO.l to 0.2. This behavior suggests qualitatively the "Reynolds number 
invariance" observed for laboratory turbulent shear layers (Cantwell 1981) . 

These slopes of 0.1 to 0.2 are even quantitatively similar to those observed 
in laboratory cases (e.g. Townsend 1976, p,181 and Table 6.5). Moreover, 

t 

the development of the large-scale non-linear structure by vortex pairing 
leads to orbit velocities ccmtparable to the eddy speeds, &v. Because this 
process is ultimately associated with the jet entrainment or spreading rate, 

'V cylindrical coordinates 2 as in CH) it is reasonable to set 

Av •• v^ s v„dR/d2 (3) 

er 2 

where v^ is the mean axial velocity and v^ is defined to be the mean 

spreading rate at the boundary. In fact R/Av ■ (8/v ) dVdR gives a good 

z 

estimate of the period of the bursting phenomenon observed in boundary layer 
flow with dZ'dR between 6 and 10, and such bursting phenomena are associated 
with the growth of the large scale structure in the layer (Cantwell 1981, p.483). 
Moreover, apparent oscillations of the radio jet axis might be expected on 

this basis with wavelength v, R/Av « R d!^dR. There are preliminary 

s 

indications (Willis et al . 1981a, b) that such oscillations exist in NGC 315 

and NGC 6251 on about the scales expected, if the jets in these sources do not 

lie too close to the line of sight. 

In a compressible jet (with the magnetic field unimportant in the 

large eddies) one expects the additional restriction Av S c , the internal 

s 

sound speed, because otherwise dissipative shocks can arise. Indeed if a jet 
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had conatant dR/dZ ovar a larga gang# of Z thia would probably indioafca an 
eaaenfcially fr«« jafe with a conatant intsmir^ found ap«td> moat plauaibly due w 
ahock diaaipation of aupersonic turbulence (a.g. BicHnell and Molroae 1981)* 

But ehould the jet be confined by the aurrounding medium (e.g. CM; BCH) then 
the turbulence will be inhibited,* and in fact the entrainment velocity 
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H Will d«cllnt to ziro with dR/dZ («■ theft ftt in ditvet proportion 
with V 2 itynptotically conatant). The tntrgy available to the intamal 
caacadt will thus alao dtclint to zero. Hanct wt prefer to work directly in 
terma of the obaerved quantity dR/dZ aa being the beat indicator of the 
phyaica of the beam. We may deduce ^ obaervationa warrant 

it (long ateady expanaion), but we may alao detect the decline and growth 
of the internal turbulence with dR/dZ '4iich provea to be eaaential when 
fitting the obaervationt (Section 1X1). 

We recall that the Kolmogorov velocity apectrum ia (i refera to the 
acale at the upper end of the inertial range) 

. ( 4 ) 


and that the effective internal viscosity of the Jet is determined when the 
Reynolds number is fixed as 

^eff (R/Ra) . (5) 

Moreover the Taylor microscale for homogeneous turbulence is given by (e.g. 
Bradshaw 1978) 

Aj “V R /I57Re , (6a) 

1 A 

from which follows 'v /I? Re , (6b) 

on recalling (2). This last equation also allows us to write ''' 

- 1/3 

(B(-5)R(2l)/(Yg(5)fh) J , which ratio is not very sensitive to our various 
uncertainties. 

We adopt the hydrodynamic view (e.g. Falco, 1977; see also Cantwell 1981) 
that the Kolmogorov equilibrium cascade exists between A^ and Aj^ as our 
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flrit; approxlnuifclon to the complex Interactloni expected In the regime, ao 
that we may set - 4^ and • Av In (4) . It is this range of eddy 
scales that we expect to Interact with the particle spectrum by way of Light- 
hill radiation. Depending principally on the effective value of fp (see 
below) this range is always 10** in scale, and thus also in particle^ energy. 

The Lighthill formulation of radiating eddies as developed for the MHD 
case by Kulsrud (1955), Parker (1964), and especially Kato (1968) is 
summarised by Stein (1981). In the presence of a magnetic field, the Al£v6n 
waves (A) and the slow magnetosonic waves (S) are emitted predominantly in a 
monopole mode. The power emitted in this mode from an eddy of scale 4 is 

■■a.s ■ I* • 

where X is the wavelength of the radiated wave, and we have written the 
emission efficiency as n • 4/X . The fast magnetosonic waves (F) are emitted 
either In a dipole mode (confined jet) or in a quadrupole mode (free jet) with 
the power (n 1, dipole; n « 2, quadrupole) 

Pp A. rij. • (4/X)^=^"'*'^^ p • (Av^ 3/4 , n - 1.2 . (8) 

In these expressions we should take (Av^)3/4 » (Av)3/r as this last 
expression is thu source of the available energy. (The time avera^ie is over 
times long compared to the Taylor ’’turnover time”. It would be incorrect to 
use (Av) 3/4^ because this power can only be maintained sporadically in ’’bursts” 
of duration (4^/R)(R/Av) - i^/Av).) 

The principal factor in the emission efficiency is determined from 
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it/\ i f • 


(9) 


where c Is the large scale AlfvAn speed, , for A or P waves, and 
the mean sound speed, c^ , for S waves. The factor f , introduced 
in equation (1), is identified here as the eddy scale in units of the wave- 
length of its principal radiated wave. 

Equations (4) and (9) allow us to write in fact that 

(Av/c)^"^^ s f , (10) 

Moreover we can estimate Av/c^ from the results of CH, Bridle al, (1980) 

■.2 _2 

and BCH, which show that the quantity Cg 'Vi (c^/V^og) litust be K 10 in 
order to prevent large-amplitude pinching oscillations (largely unobserved) 
in the regime where B is predominantly circumferential. Hence using (3) 
and c^/vg ^ 10“^ obtain Av/c^ ^ 10(dR/d7,) - 0(1) . Thus %l\ 1 on 

the Taylor microscale, and it declines like or (see (4) and 

(10) and recall (9)) towards the smaller scales. We have already argued that 
Av/c^ 1 , so that we may conclude that the monopole emission of A and S 
waves will generally dominate the F wave emission (jj, < Aiy) • i'®® also 

that the emission efficiency g ^-/X decreases as ~ *^A S 

[f(Jl)]V3^by (4) . From this expression for f(Jl) and the expression for 
we may estimate f 'V lO”^'^® tyYg(5)/(B(-5)R(21)]^/® 10“®^® . 

Consequently on inserting this value in equations (1) , (2) and (6b) we find 
-V 1013 (,iYg(5)/B(-3))cn , Re 'v- 10lO*8 [B(-5)R(21)/(yYg(5)) , and 

103.‘*[B(-5)R(21)/(uYe(5)]3/8 . It should be remembered that B R 

'V const, in these formulae (e.g. CH, BCH), so that the numerical values are 
reasonably well known (Aj^,Aj will scale with R) . 
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The decline of the tadiiitlon efficiency of a turbulent eddy with C(Z) 

Is offset somewhat, at least for the Alfven waves (Kato, 1968), by the tendency 
for to rise with the approach of the magnetic field to equlpartltlon 
(from 3.05 for a weak field to 15,5 for an equlpartltlon magnetic field), 
which might be expected on the smaller turbulent scales (DeYoung, 1980), 


Thus we mlgnt expect s^irong Llghthlll emission of MHD waves from eddies over 
the article Inertial range from to . We know moreov-r from the work 

of Lacombe (197'^') and Ellek (1979, 1981) that such waves can be damped 
locally (10^-10*’X) by acceleration of supra-thermal electrons (we suppose 
sufficiently low thermal densities that heat losses do not dominate, as this 
gives reasonable estimates for the sources studied here. In general the 
"thermal efficiency" must be calculated from the various damping rates.) 

We clearly are confronted at this point with a difficult and interesting 
calculation. We wish to know the self-consistent particle, MHD wave and 
indeed (allowing for the back-reaction) turbulent eddy energy spectra that 
are In equilibrium over the inertial range. We defer the proper solution of 
this problem to a subsequent investigation, but a useful estimate of the 
particle, spectrum can he made by neglecting the back reaction on the eddies, 
as follows. 

From Equations (7) (time averaged), (9), and (10) we may compute for a 
unit volume the flux density of the "Lighthill radiated" MHD waves, 
w(k) , as (k = 2ir/X) 


woo = . p- (11) 

where k^ = , Using the expression for the resonant wavelength in terms 
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# 


of the particle energy, namely k - l/(yrgg) - (eB/y)E”\ we also have 


dP 

dE 



( 12 ) 


We equ'^te thla damping power per unit particle energy to the synchrotron power 
dPs/dE radiated by the system per unit particle energy in accord with the 
view that this radiation is the principal dissipative mechanism. We have then 



(13) 


where (dN/dE) is the number of relativistic electrons per unit energy per unit 
volume with energy E, so that the particle spectrum should be of the form: 


dE 


CC 


^- 2.5 

E 


(14) 


for constant local conditions. This corresponds to a synchrotron emissivity 
e oc which is only a little steeper than the v spectrvan observed 

V 

for typical radio jets (Bridle 1981; Willis 1981) • When the numerical 

coefficient in (14) is computed (by equating (12) and (13)) and the minimum 

particle energy is taken to the 10^ , we deduce that the density of 

relativistic electrons n 10"®n near the bases of the jets. We empha- 

er e 

size however that, although these estimates are plausible, the assumption of 
balanced input and output over the whole energy spectrum is arbitrary and 
none of the above can be considered proven until a careful calculation of 
the various spectra (Including heat losses) has been completed. 

our estimated spectrum is approached at higher frequencies in the 

bases of the jets in 3C 66B, 3G 31 and NGC 315 however, for which Butcher et al . 
(1980) determined "radio to blue" spectral indices of 0.70, 0.72 and >0.74 


15 



respectively. This is to be expected because particles radiating at optical 

frequencies will have shorter synchrotron lifetimes and are therefore more 

likely to be in equilibrium with the turbulent input. In fact we might expect 

the spectral break to be at the particle energy for which the synchrotron 

lifetime t (y ) is about equal to the turbulent "turnover time" or "convective 
s e 

time" T 'V z/v "V r/Av. This is because a family of electrons near this energy 
c z 

can radiate energy at just the rate at which it is both convected/r and delivered 

13.7 “I “2 

to the viscous subrange. From the relations 10 * Yg <5) B (5) , and 

T 10^^*^ (z(kpc)/v (8))» we see that Yg(5) 'v 1 (mm or far IR) would 

be so selected near the base of a Jet. However, unless v * l/R , this energy 

z 

would increase as R after the magnetic field transition (bI^ “ r 2 and Z « R) . 

We also observe that if collisionless shocks are to appear in our 
scheme then we would expect them to do so near the Taylor scale 5.^ , because 
there the Lighthill waves are non-linear. These might in turn lead to a 
fraction of the turbulent energy being diverted to shock-wave Fermi accelera- 
s ion from direct resonant wave acceleration (e.g, Bicknell and Melrose 1981; 
Blandford and Ostriker 1978; Bell 1978a, b), but we are unable to estimate 
what effect this might have on the MHD wave or particle spectrum at present, 
and so we shall not consider it further here. 

In summary, we propose that the turbulent energy available from the 

large scale structure in the jet "mixing layer" is dissipated through the 
acceleration of relativistic particles or scales from to (much 

larger than the proton gyro-radius and Debye length). The detailed inter- 
action of turbulence, MHD waves and thermal and non-thermal particles in 
this regime is left undlscussed, but a particle spectrum with a 

^i®f^®uing below E/ro^c^ 'U 10®R(21) seems possible. We proceed now to 
deduce a singular observational consequence of this scheme, and then to test it. 
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(b) The Scaling Law for the Radio Surface Brlghtneag 

Given the suggeetion that the synchrotron emission is the basic 
dissipation mechanism for the turbulence in the radio jets, we can write 
immediately that 

e^/P ''' e(v)(Av)®/R (15) 


where the "jet radiation efficiency" e(V) contains the prescription for 
the spectral distribution of the energy. In general there will be an addi- 
tional factor for thermal efficiency say (1-e^) , where Is the frac- 

tion of wave flux density w(k) lost to heat, and may be quite close to unity 
in some cases, Hc'ever, we pursue our argument here with e^«0 for clarity. 

It is reassuring ph* ;ically that an expression of the form (15) may be calcu- 
lated directly from equation (12), (13) and (14) (with the various coefficients 
collected and using the standard formula for the synchrotron emissivity of a 
power law: e(v) a V * ) , but in view of the various uncertainties in that 

procedure we prefer to base this section on simple scaling relations. Our 
basic idea does require (and the estimates of the preceding section justify) 

taking (neglecting losses to the thermal particles completely) 

V 

I* max 

dv - l/4ir , (16) 

o 


and, if e (v) oc v 


-ot 


then 


^ -J* ®('^) 


4tt max 


(17) 


- 0 . 6 + 0. 2 

As the observed spectra of most radio jets are “ V (Bridle 

1981; Willis 1981) we take the jet to be optically thin, so that the surface 
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brightness B(v) c(v) • a , where s is the thickness of the jet along 

the line of sight (/r2"Iso'^‘ at a perpendicular distance -'fij from the centre 
of a circular Jet). This gives a geometric limb-darkenlog for the jet such 
that, in the absence of intrinsic emlssivity variations, the half-power 
points are on the circle 'dS > (✓3?2)R , Such limb-darkening is consistent 
with the observed transverse intensity profiles in 3C31 and NGC315 (Fomalont 
et al. 1980). We now combine equations (3) and (IS) with an equation for the 
mass flux, A , in the form (we consider a variable mas(i flux due to entrain- 
ment in Section IV). 


R^p V ■ A 
z 


(18) 


-2 -1 -1 -1 

to write the surface brightness (erg. cm s sr ) as 


B(V) * (1-(S/R)2)1/2 e(v)A (v /R)2(dR/dZ)3 , (a) 

z 

that is (19) 

B(V) (l-Cffi/R)2)l/2 e(v)A v^/(r2v^) . (b) 


In practice/, B(V) is the quantity determined by observations of a 
well-resolved jet at a single frequency v • However, although the form 
of the single frequency brightness variation should be given adequately by 
(19) the amplitude is more problematical because of the uncertainty in e(v) . 
We feel that the Integrated result (16) is more reliable in this regard so 
that we will use the integrated radio brightness for amplitude determinations, 
namely 

B = (1-('5)/R) 2)^/2 e A (v /R) 2(dR/dZ) ^ , (20) 

z 


where 


« 1/4TT (V_/V )<^““^ 


CO max 
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by (17). With the radio band "cut-off" at 'v 10^® Ha, and 4xl0^' 

(Yg -6 lO^Aire » 0,4 for a • 0.75, and M).l6 for <x 0.5 . Even if the 
synchrotron spectrum extended Into the X-ray band K 10®) , we would still 
have of the total turbulent power In the radio band for 0 0.75 . 

Wlilch of these various energies actually obtain probably depends on the energy 
of Injection of the particles at the base of the jet, but at least in the low 
luminosity sources the Integrated radio brightness Is likely to be a sensitive 
Indication of the total turbulent power. 


m, COMPARISON WITH OBSERVATIONS 

In this section we compare the predictions of our basic scaling relations 

given in equations (19a) and (19b) with the observations of two well-resolved 

radio Jets, those in NGC315 and 3C31, The variations of jet radius R with 

distance Z from the radio cores have already been modeled for these systems 

by Bridle etjadj^ (1980) and by BCH; these models provide us with estimates of 

the variation of v with Z , essentially by requiring that the flows satisfy 

z 

Bernoulli's equations 
(a) NGC315 

A comprehensive study of NGC315 at 610 MHz and 1415 MHz has recently been 
made by Willis et al. (1981a, henceforth WSBF) . Their Figure llA shows the 
behaviour of the l4l5-MHz surface brightness B(v) on the axis of the jet 
as a function of its PWHM. These data are replotted in *lgure 1. The curve 
superimposed on the data shows a prediction from our equation (19a) based on 
the BCH fit to the collimatlon variations of the jet. This procedure essen- 
tially adds the form of the variation of v^ , which is known from the BCH 
fit, to the smoothed values of dR/dZ and R(Z) observed by WSBF. The 
brightness normalization can be considered arbitrary for Che moment. 

There is clearly excellent agreement between the shape of the curve pre- 
dicted from equation (19a) and the observed brightness variation. This is 
encouraging, as in the absence of field or particle replenishment in the jet 
the surface brightness would fall off as R , as shown in Figure 1 by the 
dotted line. Furthermore, on this logarithmic scale the fit is not very sensi- 
tive to details of the v^^ variation determined from the BCH fit. As discussed 
by BCH, a variety of models with different combinations of pressure and magnetic 

% 

confinement could be fitted to the same R(Z) data to within the observa- 
tional errors; the goodness of fit in our Figure 1 is insensitive to this 
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r«nge of alternativei. 


We note that equation (19a) succaaafully eimulataa the Initial regime 
in NGC315 wherein B(v) falls off as r“^* 25±0.12 according to WSBF. 

Inspection of equation (I9a) shows that when dR/dZ is approximately con- 
stant, and there is no significant entrainment, we expect the brightness to 
vary as v|r”^ , so that the slowest Intensity decline In an unaccelerated 
regime would be r“^ . The ability of the model to reproduce the data 

therefore stems In part from the fact that the BCH fit to the colllmatlon 
data required this to be an accelerated regime in which v Increased by a 
factor of about 2.5 . H»e success of the fit to this regime is therefore 
intimately connected with its proximity to the ’’sonic height” of 5” (1.2 kpc, 

Hq - 100 km s"^ Mpc'^b in the BCH model for the collimation. Turbulent 
entrainment would provide an alternative explanation (see below), but the 
collimation data suggest a rapidly declining ambient density which would 
inhibit entrainment asymptotically. However, small variations in A can 
occur near the base even then (see (25) in Section IV below). 

This Initial regime is in marked contrast to the regime beyond 2R 'V 20” 

and extending to 2R 'v 45 *' . Here there is little variation in v^ in the 

BCH models and the goodness of fit displayed in Figure 1 reflects the smoothed 

. 42 ± 1 .10 

approach of dR/dZ to zero at 2R 'v 45” (yielding B(v) « R •• (WSBF) - 

hence R“ Z®/® on this shoulder by (19a)). Subsequently the predicted 
turbulent power drops to zero, but the observed brightness is sensibly 

instant over this whole flat region between 2 *.5 and 7 ’ from the nucleus 
(where 2R Oi 45 ”) . 

In the outermost regime the predicted turbulent power rises smoothly 
with dR/ dZ to the correct observed level above the adiabat extrapolated 
^^*^^^ateau and then declines in step with the observations (approxl- 
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mately at r“* (WSBF), Once again thete it no need to Invoke efficient 
entrainment in thit outer expantion, and the dentity law obtained from the 
coUimation fit (BCH) would not predict it by (25). below. Thut over the 
whole well-obaerved main let in NOC315. the predicted briahtneta tracks the 
obacrvatlona remarkably well , excepting only the collimation plateau. This 
hat certain Implications for the physical parameters of the jet as we now 
proceed to demonstrate. 

Firstly* in order for the synchrotron eroisslvlty to vary in proportion 
to the local turbulent power input as in Figure 1, the Intermediate resonant 
waves must be strongly damped locally, pilmarlly by acceleration of relativis- 
tic electrons. That is, the presumed acceleration must be both efficient and 
rapid (as a convective time scale). This requirenmnt is in accord with the 
damping scales estimated by Eilek (1981) and with the dominant damping mode 
expected (supratherraal acceleration rather than heating; e.g. Lacombe, 1977). 

Secondly, ’’convective smearing” of the turbulence dominated B(R) rela- 
tion will be avoided only if the lifetime of the radiating particles la 
shorter than the local convective or ’’turnover” time scale (t 'V AZ/v a 

C 2 

R/Av) . On the collimation p‘’atea.} the synchrotron lifetime of the electrons 
radiating at 1415 MHz „n the equipartition field strength of 4pG (WSBF 
rescaled to Hq ■ 100 km s Mpc ) and in the presence of the 2.7K background 
would be Tg "V 7.6x10^ years. The length of the shoulder region is Az 'v 65 
kpc so that V need only be greater than AZ/t„ “v 900 km s"'^ for convective 
’’smearing” to occur. This is quite likely according to our estimates below, 
and so we expect the observed profile to be flat in this region, as is observed. 

In regions such as the inner and outer expansion regimes where the bright- 
ness predicted by equation (19) is above that of the expansion adiabat, 
the local turbulent input will dominate 



as obsQ^vsd (figuv® I) ptovidad only that it it dissipatttd ot oonvairt«d on a 

conviction tiroa seal® ot loss. Tha integtated (10 mHjs to 10 GHz) xadio 

“ 0.6 

bxightniSSi assuming that tho spsctxvun is v over this range 

(WSBF) should than to® given toy (20) with 4ff® « v”®j^^{lQ) i 1 /h» We now 

determine the scaling coefficient in equation (20) toy fitting it to the 

absolute brightness of the jet in NGC 31S. 

we equate the observed and predicted brightnesses at R « 2” (472 pc# 

■*1 *•! 

Hq » 100 km s Mpe ) where the jet ia well resolved toy the VUV. using 

A R 

equation (20) with *5T ■> 0 then requires that B/eA » 8.9 x lo"'’" v| where 

“1 

v^ is the flow velocity at the "notsla" in an s . The observed radio 

brightness at this part of the jot (interpolated from table 3 of llSBF, 

“6 “2 “1 -1 

renorm-aliced to « 100) is 5.2 x 10 erg cm s sr , so that we require 
38 

eAv*^ n» 5.9 X 10 e.g.a. 
s 

Wt^ now use 4ita s \iax^ accord with our earlier argument (the 

inequality reflecting our uncertainty about the energy lost to "heat")/ 
together with equation (18) / to obtain directly (and relatively firmly) 
that! 

p Vg v| ^ 3.5 V°;^(10) X lO**^ e.g.a. (2J) 

where p and v are the density and longitudinal velocity of the flow at 

the distance from the nucleus where the jet radius is R *» 472 pc. These 

parameters can be rei'laccd by parameters at the "nosale" using the DCH fit 

to the collimation behaviour of the jet, which gives v^ » 2.1 and 

P « Pg/15* Consequently, recalling also that v| « YPg/Pg a>'d taking 

Y « 4/3, relation (28) can be written in terras of the pressure p and jet 

velocity V at the "noasle''; 
s 

-2 0*4 . -2 «1 

p V t 1.9 X 10 V (10) erg cm s 
S s max 
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One ieee Inmdlecely from thii eitlmate of the energy flux in the Jet together 

with the BCH value for the Jet radiui at the norrle ■ ISOpe) that the 

lob. lu.lno.lt). L n, 4 2 (10> 10-» .r>. . Thl. v.lu. 

L v”®***) le iatlefaceerily cloae to the obaerved radio lumlnoaity of each of the 
miix 

north-preceding and aouth following lohea (Bridle et 1976) » hut thta 
cannot be regarded too aerloualy In view of the varioua uncertaintiea (expec- 
ially in v^,, R^ and e ) . 

We may alao use the BCH fit to acale p^ from the equipartition calcula- 
tion of WSBF. Thii givea (aee e.g, Parley « *1*. 1979 for brishuneaa depen- 
dence) p^ A. 1.5x10“^ (Ih) dynea cm"^ ao that 

Vg ^ 1300 km a"^ . 

and (22) 

Pg l.lxlO"^® gm cm"^ 

With the BCH value for R (180pc) theae valuea yield a mass flux of 

s 

2x10 ^ ^ which is reasonable. 

Moreover I on the collimation plateau | (where R 'v 22" (5.2 kpc)and 

dR/dZ A/ 0), WSBF concluded from their depolarization data that the thermal 

«>4 »3 

particle density must be ^ 4x10 cm . Using equation (18), the above 

value for the mass flux (xl/t) , the BCH value for v^ on the plateau 

(“ 2.55 v^ ^ 3400 km a"' by (27)), and the observed radius we calculate 
s 

n 3xl0”®v“®’^** (10) cm”^ . Such a low value is compatible with the 
e max 

measurements of WSBF in that their data were Just consistent with no depolari- 
zation between 610 MHz and 1415 MHz. Thus the mass flu? : raguired to deliver 
the turbulent energy necessary for particle reacceleration in NGC315 is well 
within available constraints on the thermal mass in the observed jet . Should 
a thermal density much larger than our estimate be confirmed, entrainment 
occurring just before the collimation plateau is the most likely explanation 
(see Section IV). 
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For It ii thirt that th« coUiimting prtiiurt profili (and charefora pottlbly 
tiha daniicy proflla) ii fiaccait according to tha CH fit glvan in BCH. 

(b) 3C31 S N6C383 

Ttia dlitrlbutlon of Intanalty with diatance from tha radio core in tha 
Northern Jet of 3031 haa bean obaarved at 4385 MHk and at 1480 MHz with 
the (Incomplete) VLAby Fowalont ct al» (1980) and there are aavtral pubXiahed 
CH flta to tha collimation behavipur (Bridle at al, 1980) , Unfortunately the 
jjeti in thia aource undergo large deflectlona ending in the lobe atructure 
before reaching their firat regime of dR/dZ« 0 * Conaequentlyi fita to the 
collimation propettiea uaing the CH/BCH formulation are not so tightly con- 
strained as they are for NGC315, 

Figure 1 shows the variation of surface brightness, observed in the 
Nor the nt jet of 3031 at 4885 MHz , as a function of the FWHM of the jet. 

This plot is obtained by condjining the data of Fomalont et al, (1980) with 
that of Bridle et al. (1980) . Curves A and B are the predicted brightness 
variations using the collimation model (a) of Bridle et al. (1980) for A, 
and a pressure dominated OH fit (whose parameters are given in the figure 
caption) for B. 

The most important variable affecting tmc goodness of fit to the observed 
brightness variation is the height of the sonic point . Tlie best agreement 

between the predictions of equation (19) and the data for 3031 is obtained 
with sonic heights in the range I”. 5 < Z < 2'*, 5 (0.35 to 0.67 kpc) . This 
is comparable to the height Inferred for the nozzle in NGC315 by BCH (1.2 
kpc) and is for both sources remarkably high* Figure 2 also shows, however, 
that the predicted B(R) variations for 3C31 begin to decline rapidly much 
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•oontr th«n cht d«t8| lo tnat iht pradictid turf act brighcnatt In tht nora 
axpanded ragimat of cha 3C3I jac fallt balow fhae obtarvad by a factor 'VIO. 
Navarthalatti it it clear that our tnodal la uuch clotar to the obaarvad 
behaviour In the inner jat saglma than la the adiabatic axpanaion law (ahown 
dottad in tha figure). 

iha r<?giaa of rapid dacline in tha predicted brightnaaa relative to tha 
obaarvad brightnaaa (and indeed to the adlabat) bagina at 2 <v 20*' (correa- 
ponding to tha FWHM 5"). Juat aa In NGC315, this la preciaaly where dR/dZ 
bagina itta decline (the point of inflexion d^R/d2^ « 0 } tea a.g. CH) to 
what would probably be a colllmotion plateau if the jet remained undiaturbed 
(it may liappan beeauee the jet la diaturbed at thia dlatance if the diaturb- 
Ing effect Impliea a preaaure plateau; aee below), Thua we have in effect 
for 3C31 only the firat segment of the two-piece curve shown for NCG315 in 
Figure 1, and it is dropping below the data just aa does the initial segment 
of the NCG315 profile. 

The redahifta of NGC315 (.0167) and of 3C31 (.0169) are almost identical 
so that we may compare the angular scales directly. On this basis the convec- 
tive timescale for 3C31 on the shoulder of the plateau is about 1/5 that of 
NGC31S (at the same jet speed), but the equipartltion magnetic field is some 
ten times greater (MxiO’*® gauss) so that the synchrotron lifetime is 100 
times smaller at the same energy (there is another multiplication factor 
•vl,? in t to allow for the difference in frequency). Therefore, unless 
the jet velocity in 3C31 is some 12x900 km/s (sec the argument for NGC315), 
or unless the equipartltion fields are inappropriate in this region, we should 
expect convective smearing of the turbulent power to be reduced in this source. 
In fact the observed B(R) variation (Fomalont et al, 1980) continues rela- 



tiveXy flat, A relatively high jet velocity of >10** km e“* might not be 
Implauaible (aee the eetlmatea made by Parley et al, 1979 for 3C 449). 
Moreover, if there were a density plateau, and not merely a pressure plateau, 
associated with the slowing down of the lateral expansion of this jet, entrain- 
ment (Section IV below) could prevent the turbulent power from dropping so 

rapidly with dK/dZ , and thus lead to better agreement with the observations, 
Th^ 4885 MHz brightness at the base of the jet in 3C31 has been measured 

in to about 2" from the radio core (Fomalont et al, 1980) . Their Figure 3a 

shows that the jet brightness initially rises to a broad maximum about 4" 

from the core '.^fore the subsequent decline along the large-scale jet. This 

initial rise is fitted approximately by curves such as A in our Figure 2 

when these are computed to smaller values of Z using the CH models for 

the jet collimatlon. (The actual variation of jet radius with Z in this 

regime is difficult to deduce as the jet is barely resolved in the transverse 

direction by the observations of Fomalont al. 1980) . It is significant 

that this same initial rise and fall in brightness of the jet is also seen 

at optical frequencies (Butcher et al. 1980; van Breugel 1980), implying a 
-. 0,72 

V "radio to blue" spectrum. These observations in the optical band 

(which cannot be subject to convective smearing) follow the turbulent power 
through a predicted rise and fall immediately beyond the "nozzle". This is 
encouraging evidence for the basic principles of our model. 
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IV. ENTRAINMENT AND THE MASS FLUX, A 


We have seen above that, principally because of the probable huge Rey- 
nolds numbers in the observed jets, we expect tliat the observed radio jets 
are turbulent Jets . Moreover, this hypothesis leads directly to equation (19) 
which proves to be rather successful at fitting the B(R) variation of at least 
some radio jets (§III) with the mass flux A taken constant. But turbu- 
lent jets grow by continually entraining surrounding material into the tur- 
bulent region, effectively at the velocity of the largest eddies (see equation 
(3) and following; also Cantwell 1981), and so to generalize equation (19) we 
must estimate the possible form of a variable mass flux A(R) . 

We adopt a self-similar model of a compressible, turbulent jet in a 
non-uniform background. For simplicity we may suppose the jet to be non- 
rotating about its axis and to be non-magnetlc, although neither assumption 
is essential to the form of the solution. We model the turbulent jet by 
using the Navler-Stokes equations for the time-averaged mean quantities and 
a kinematic viscosity of the form 

V ■ V r sinG , (23) 

where v is a constant speed and r,6 are spherical coordinates having the 
jet axis as the polar axis. This simply states that the dominant viscosity 
is the "eddy viscosity" in a turbulent jet. An excellent estimate for the 
speed V is furnished by (3), namely v ■ v^ dR/dZ , where v^ H v^ (9“0 
and R(Z) is, as usual, the radius of the jet at any section (« Z for 
turbulent jet, Landau and Lifshitz 1959). The jet is supposed to be steady 
in the mean (t > = r/v^) . 
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This model yleXds Che following form of solution to the equations of 
motion (Henrlksen and Wang 1981) 

p - Pg r“" p(0) , r 5 r/r^ , 

p - py2 P<e) , ( 14 ) 

- V v^(6) , Vg - V Vg(9) , 

where r , p are evaluated at a standard reference point; p(0) , P(0) , 
Vr(0) • dimensionless functions to be found from the equations of 

motion (2) , the continuity equation, and an energy equation (note that 
kT(0)/(p^nijjV^) i P(9)) . We are Interested here only in the continuity 
equation which yields (R = r sinOj , 9^ Is jet half -angle) 

/ \ 

A = R^p v^ « const. - (sin0^)^” ‘ (f^ nf 

or, for n«2 ^25) 

■ const. - (R^ p v) A(0) iln R . 
s s 

where A(0) = ^ <P(9) Vg(9) sin0)]g_p < 0 . 

This last result gives a reasonable estimate of the forms A(R) to be used 
in (19) , for an arbitrary ambient density power law (n) . We may conclude 
that core entrainment will cause A in equation (19) to increase signifi- 
cantly with R if n 2 (A “ R^“” asymptotically with R). Otherwise, 
we see that A will soon tend to a constant with R . This result is 

Important, as it implies that a jet brightness that declines less rapidly 

_a 

than R (see 19a when v^ , dR/dZ are constant) may be associated with a 
radially flat ambient density profile (B(v)® R” ^ by (19) and (25) if n-1 


B(\i) 'y const, if p 'v const, (i.e, n»0 )). Should this happen near the base 


of the jet, this mechanism would explain flat B(R) variations without invok- 
ing a nozzle to produce an Increasing • The height usually attributed to 
the "nozzle" (CH, BCH) might in such cases simply be the onset of the "flat" 
density variation. Moreovera a jet which encountered density "clouds" would 
show considerable brightening or a "hot spot". It is also likely that vor- 
ticlty and density irregularities will be entrained by the jet from ambient 
density "clouds" (also including the nozzle region) and convected outwards. 
This will introduce a "hot spot" structure in the jet. For a sufficiently 
extensive density plateau one might expect a periodic hot spot structure on 
an Interval RdZ/dR . 

We observe however that, for a given beam energy, the turbulent beam 
must develop a core-halo structure (Vj.(0) , Vg(6) , and possibly p(0) 
must be peaked sharply near 0*0 with broad "wings") to compensate for 
the additional mass being added to the beam. The numerical demonstration of 
these conclusions must be left to a subsequent paper, but we note that our 
remarks here are generally in accord with the numerical results of Baan 
(1980), who used a similar turbulent viscosity. 

Finally, in this Section we give the .analogues of equations (24) and 
(25) when there is an effective viscosity, V , which is constant. These 


are 


p - pg r"" p(0) , r s r/r^ , 
-(n+2) 

p - Pg -pr P<9) ^ • 

s 


V 


V ''0 


VflO) 


r r 


'0 r r 
s 


(26) 
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and 


and 


A “ const. - 
n - X , 

A “ const. - 


vp^ R^(slnep” 

(X-n) 


^^“"^A(0) 


(vp R ) 8in0. (in R)A(O) 

SB 3 


( 27 ) 


We see that the core entrainment is Xess rapid for a given density profiXe 
than in the turhuXent case (but the totaX entrainment may be greater due to 
a fXatter v^(d) depending on the nature of V : see aXso Baan X980), but 
that the jet is effectiveXy deceXerated. This might be the fate of inltiaXXy 
reXativistic nucXear jets untlX they become turbuXent. but the nature of the 
necessary viscosity is probXematicaX (Baan X980) . 
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V. CASCADING JETS 


The nature of the connection between the nuclear (VLB) jets and the 
halo (VLA) jets must be addressed in view of the sources observed to have 
both. Early opinion (see e.g. Miley 1980) seems to have considered the large 
scale Jet to be a simple extrapolation of the nuclear Jet, but this leads to 
enormous nuclear sources of magnetic flux and possibly angular momentum (if 
the CH helical fields are supported). The resolution of this problem would 
seem to be the entrainment of mass, magnetic flux and angular momentum by the 
beam as an inevitable part of the cascade to the larger scale (CH and De 
Young, 1980). One must, however, reconcile the similar position angles of 
the VLB and VIA Jets that are observed in at least some cases (e.g. NGC315 
and NGC6251), with the mechanism of this cascade. We attempt in this section 
to summarize the various possibilities. 

Should the ambient pressure, p^ , in the clrcumnuclear material be such 
that p^ « pj vj^ (pj and v^^ refer to emergent properties of the nuclear 
jot), then we would expect the jet to be relatively unrefracted (Henrlksen 
e^t al. 1981) . When p^ Vj^ deviations from the initial direction 

can occur however, unless this initial direction is accurately aligned with a 
symmetry axis of p^ . The case p^ > p^ v^^ will not arise in an existing 
jet, because it would produce Instead a confined core source. 

In order for a high "nozzle” to form by the Blandford-Rees mechanism 
(CH, BCH and this paper), and for a memory of the "original" (after nuclear 
refraction as above) direction to be simultaneously preserved (as suggested in 
CH), the material entrained by the nuclear jet in any given interval must also 
be significantly heated. For only then can a "bubble" of heated gas, which 
retains a memory of the inltla' momentum, but which is also subsonic, be 
produced on the new spatial (and mass) scale. 
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This clarly requires for a transverse spatial scale, £> , that « )l^/a , 

A 

where a is the effecrive thermal diffuslvity (we always assume that the jet 
is not radiating most of its energy), Hence, the Prandtl number ot the nuclear 
entrainment = v/a should be e i . if this ratio is much smaller than 
unity then the expanded jet will not retain Its original momentum sense, while 
if the ratio is much greater than unity the jet will not renozzle. We must 
also suppose that the cooling time of the heated gas is larger than , in 

order than a "fire box” for the nozzle be in fact established. 

Given that P^ » I and a sufficiently long cooling time, whether signifi- 
cant entrainment does in fact occur may be estimated in general from equations 
(25). However, a condition which is always sufficient (but usually not neces- 
sary) to ensure the jet*s growth is simply JL*‘/y « z/v^ , or Re = * z/£ 

^ cot 0. . 
j 

We conclude therefore that v “ a 4 sufficient conditions for an 

"upward jet cascade with memory”. Moreover, the easiest wav of satisfying these 
conditions is by a turbulent transfer of heat and momentum from the jet to the 
ambient medium . We are led then to a model in wMch the nuclear jet disperses 
its energy and momentum over a much larger mass scale by inducing vortical 
turbulence in the ambient nuclear atmosphere which results in efficient heat 
transport and entrainment. This nuclear "fire box” region can than renozzle 
to produce a jet of lower specific energy but carrying in exchange much more 
thermal material as well as angular momentum and magnetic field (the angular 
momentum is likely to be present in the ambient material, while the magnetic 
field can be produced in part by the induced turbtslence, DeYoung, (1980)). This 
"fire box” region may be X-ray bright as a thermal source but, without the 
energy input from a mean motion, as discussed in Section II, we could not 
expect a strong inertial cascade to be maintained nor thus a strong synchrotron 
spectrum. 
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VI. DISCUSSION AND CONCLUSIONS 


We have attempted in Section ll of this paper to render plausible on 
physical grounds the proposal that the synchrotron radiation of relativistic 
electrons provides the effective viscosity for radio jets , and that these jets 
should be viewed as turbulent mixing regions . This proposal has still to be 
justified in detail theoretically i but it is shown here to lead to remarkable 
agreement with observations of the radio brightness in NGC 315 and 3C 31. 

Our predicted brightness-radius variation (19) , although subject to 
uncertainties (which we have discussed) regarding entrainment and spectral 
distribution, is based on very general considerations, essentially on scaling 
relations. It is therefore not expected to depend very sensitively on the 
details of our theory. In Section III we have argued that the VLA and wsR'r 
data on NGC 315 and 3C 31 strongly support our prediction, as the expected 
surface brightness variation with jet radius fits the observations of NGC 315 
well, and fits that part of the observations of 3C 31 where it may best be 
expected to do so. 

Our prediction relates the local synchrotron emissivity in a jet to the 
local turbulent input , and does not apply in regions of the jet where the 
emissivity can be maintained by convection. Such regimes become apparent in 
both NGC 315 and 3C 31 where dR/dZ 0. Our models are self-consistent in 
that they do yield the convective time as less than the synchrotron lifetime 
in these very regions. 

We have also made a preliminary fit to the absolute brightness of the NGC 
315 jet, thereby determining the unknown amplitude coefficient (e(V)A) in our 
equation (19). The deduced mass flux and velocity range in the jet appear 
very plausible. We have not found it necessary to include ' turbulent entrainment 
in our fits to NGC 315 or 3C 31. This has required substantial longitudinal 
acceleration of the flow near the base of the jets to fit the slow variation of 
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brightness with jet radius there (regions of large dR/da) . Our model for the 
brightness variation thus reinforces the concept of a ''high noaale'* that was 
suggested independently (from collimation data) by CH (in their model and for 
3C 449), by BCH for NGC 315 and by Bridie et al . (1980) for 3C 31# The rapid 
decline in ambient density required to produce the Initial rapid expansions 
of these jets at their bases is also consistent with neglecting entrainment 
according to the estimates made here in Section IV. 

We cannot exclude an alternative model however, wherein the jets have 
encountered density plateaux at their bases and efficiently entrain material 
as they grow (Section IV) according to (24) and (25). This entrainment would 
continue only until tlie ambient density variation changed its form drastically, 
such as during a steep decline from the nuclear atmosphere or in the collima^ 
tion plateau (dR/dZ! ’*■ 0) region, thereby breaking the similarity. According 
to the ideas expressed in Section V (elaborating those of CH) such a hot region 
might itself be the "cascade" or "nozzle" region wherein the nuclear jet is 
growing to the scale of the "VLA" jot. 

The entrainment is also likely to be relevant to the existence of "hot 
spots" in the mean behaviour of the beam. For should the jet encounter a 
density "cloud" in the average atmosphere (of sufficient extent) which it 
entrains, then the amplitude of (19) will increase anomalously and a "hot spot" 
will form. 

In any event, both the B(R) variations and the collimation behaviour of 
NGC 315 and 3C 31 (and possibly 3C 449? CH) lead us to conclude that the "VLA" 
jets are influenced by an ambient medium at heights of 0.5 to I kpc from 

the galactic center, as though the "nozzle" envisaged by Blandford and Rees 
(1974) were there. This raises the question of the relation between the 
"VLB" jets and the "VLA" jets, particularly when they agree very nearly in 
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position angle. NGC315 itself provides a clear example of the problem, Lin- 
fleld (1981) has detected a jet about 0".003 (0>7 pc) lung extending from the 
unresolved core of NGC3A5. Its position angle agrees to within errors (and 
Intrinsic oscillations) of ±3*^ with the position angle of the base of the 
"VLA" jet. Tliis "VLB" jet rapidly fades in brightness with increasing dis- 
tance from the core, wiiereas the "VLA" jet rapidly increases in brightness at 
the very distance above the core modeled by BCK us its "sonic height 
(Bridle et al. 1979; Bridle, Fomalont, Palimaka and Henrikaen, in preparation). 
These data taken together suggest that the colllmation of the large-scale 
("VLA") jet in NGG315 has been at least a two-stage process. The first stage, 
or "trigger mechanism" occurs on scales less than 1 pc, and produces a jet 
whose rclacivistie particles are rapidly cooled » the quickly-fading "VLB" 
jet. We suggest (as did CH) that the interaction between this high specific 
energy (probably relativistic) jut and a surrounding gaseous atmosphere in the 
inner kpc of NGC315 has permitted the formation of a second "nozzle" (as 
discussed in Section V and above). Moreover we maintain that this strong 
interaction has restarted a turbulent cascade capable of accelerating elec- 
trons to relativistic energies in a regime (the "I^A" jet) where the cooling 
time for the particles is now much longer. 

We are making further observations of the brightness distributions at the 
bases of large-sr-le radio jets using the VLA to examine whether there is 
similar evidence for "second-stage" reacceleration in other radio galaxies 
on kpc scales. 

Finally we note that the model for the brightness variations presented 
in this paper must eventually be reconciled with the observed linear polariza- 
tion amplitude and structure. The amplitude has values exceedii>>; 60% in the 
outer it'*ln jet of NGC315 (WSBF) . Both this striking fact and the polarization 
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■tcuctura have been InCerprated alf.heir in terms of largn scale helical fields 
(e.g. Fomalont et al. 1960 and CH) or in terms of a tangled field which lacks 
a radial component (taing 1981). It is not immediately clear whether or not 
either of these mean field structures are compatible with the ring vortical 
turbulence proposed here. We are however somewhat hopeful|because a mean 
(entrained) angular momentum coupled to the mean jet momentum will inevitably 
tend to produce a helical mean field (see e.g. CH)^ provided only that the 
amplitude of an eddy magnetic field decreases sufficiently rapidly with its 
scale. 
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FIGURE CAPTIONS 


Figure 1. The 1415 MHz surface brightness data points <ln mJy/arcsec^) 

for NGC315 are from Willis et al , 1981a. The smooth curve Is 

the prediction of equation (19a) with A constant* /7ST ■ 0 , and 

V. taken from the fit to the colllmation data by BCH. The 
z 

scaling quantity eA is the only free parameter. The adiabatic 
extrapolation from the initial brightness is also shown for 
comparison. 


Figure 2. The 4885 MHz surface brightness data points (in mJy/arcsec^) 

for 3C31 are taken from Fomalont et al . 1980 and from Bridle 

et al . 1980* Curve A is the prediction of equation (19a) with 

A constant* > 0 , and v taken from colllmation model (a) 

z 

of Bridle et al. 1980. Curve B is similar except that v is 
found from a pure pressure colllmation model wltu (see e.g. CH* 
BCH); m“5, m'«3*R ^ 0”.21 * z " 2". 3 * z ■ 10.8 * 

H ■ 20 . The adlabat appropriate for the spectral index of this 
source is also shown. Curves A*B have only the one free parameter 
eA . 
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mJy /arc sec^ 



B(48G5 MHz) 
mjy / arc sec2 
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